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Abstract: We perform a quantitative study of neutrino phenomenology in the framework 
of minimal supergravity (mSUGRA) with grand unified theory (GUT)-scale tri-linear lep- 
ton number violation. We show that only two non-zero GUT scale lepton number violating 
parameters and three charged lepton mixing angles are sufficient to account for current 
neutrino oscillation data. This allows collider studies to be performed in a manageable pa- 
rameter space. We discuss some phenomenological consequences of the models, including 
tuning issues. 



Keywords: supersymmetry, neutrino physics 



Contents 

[ll Introduction |l| 

[2]. Neutrino masses and mixings § 

|3|. Numerical procedure 12 



g. Results [15 



|5[ Discussion 22 



|AL Mass Matrices [25 



a. One loop self energies of the neutral fermions 27 



Mass insertion approximation in CPE-CPO cancellations 28 



EL Calculation of I 1 -> / J 7 [30 



1. Introduction 

In general, the minimal supersymmetric standard model (MSSM) contains baryon and lep- 
ton number violating (LNV) operators in the superpotential [1]. This generally leads to 
fast proton decay beyond experimental limits, unless an additional symmetry is imposed 
on the theory. The most widely studied symmetry is R-parity, which forbids all problem- 
atic renormalisable operators and makes the lightest supersymmetric particle (LSP) stable, 
providing a potential dark matter candidate. However, there are still dimension five opera- 
tors which may be dangerous for proton decay [2]. If one instead imposes an anomaly- free 
proton hexality discrete abelian symmetry Pq [3], such dimension five operators are forbid- 
den. Another alternative to ensure slow enough proton decay is so-called baryon triality, 
which is an anomaly-free symmetry [4] and forbids baryon number violating terms past 
dimension five. The model allows a subset of all available MSSM R-parity violating op- 
erators [5], namely the LNV operators. LNV operators in turn lead to the generation of 
neutrino masses. Thus, the LNV MSSM provides an alternative neutrino mass generation 
mechanism to the see-saw mechanism [6,7], in which gauge singlet right-handed neutrino 
superfields are added to the Pq— conserving MSSM. Although the LSP becomes unstable 
in such schemes, it could still form a dark matter candidate if it is the gravitino, since the 
decays are slow on cosmological time scales [8]. Here, we shall neglect the dark matter 
relic density, since it could instead originate from a hidden sector. The aim of this paper 
is to find and investigate simple GUT-scale LNV MSSM models that can give rise to the 
observed neutrino mass and mixing pattern. 
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The chiral superfields of the MSSM are charged under the standard model (SM) gauge 
group G SM = SU(3) C <g> SU{2) L ® U(1) Y as 



Q 

L 
E 



(3,2, i), U ■ '■' 1 2 ' n . i i 

(1>2, — |), 



(3,1,-|), D: (3,1, 3,. 



(1,2,-i), (1.1) 
(1,2,1). 



(1,1,1), H u 
The full renormalisable LNV MSSM superpotential is given by 

W = W RPC + W LN v, (1.2) 

where Wrpc is the superpotential in the R-parity conserving (RPC) case, and Wlnv is 
the LNV superpotential given by 

Wrpc = (YE^jHaLiEj + (Yn)ijH d QiDj + (F^)„Q i iT u f7 j - fiH d H u , 

Wlnv = ^ijkLiLjEk + X' i j k LiQjDk — mLiH u , (1.3) 

where we have suppressed gauge indices and {i,j, k} £ {1,2,3} are family indices. (Ye), 
(Yd) and (Yf/) are 3x3 matrices of dimensionless Yukawa couplings; \jk, A^- fc are the 
dimensionless tri-linear LNV couplings, and [ii are dimensionful bi-linear LNV parameters. 
The soft supersymmetry (SUSY) breaking Lagrangian with LNV is given by 



C so ft = £rpS + £rpc + £lnv, (1-4) 



where 



£rpc = \MiBB + ±M 2 WW + \Mzg~9 + h.c. 

+Q ] {m^)Q + fj\ml)fj + D^(m 2 & )D + V \m\)L + &{m\)E 
h%, 

£ % RPC = {hE)ijhdLiEj + (hr))ijhdQiDj + (hu)ijQih u Uj — Bhdh u + h.c, 
£lnv = -jhijkLiLjEk + h' i: j k LiQjDk — L>iLih u 

+h\m\ iL .L i + h.c. (1.5) 

where tilde denotes a super-partner of the more familiar Standard Model field. Mi, M 2 
and M3 are the masses of the super-partner of the SM gauge bosons, the m 2 's are the 
scalar mass parameters of the sparticle and/or higgs fields, and the h's, B and Di are 
trilinear SUSY breaking parameters that correspond to the dimensionless supersymmetric 
parameters displayed in Eq. ( |1.3| ) . Tree level neutrino masses originating from such a theory 
were derived in Refs. [9, 10]. Some of the weak-scale parameters in Eqs. (1.3,1.5) were 



bounded by neutrino oscillation data in Refs. [11,12]. The bi-linear couplings have been 
employed to fit solar and atmospheric neutrino oscillations at the weak scale [13-15]. Work 
on bi-linear couplings in the MSSM can also be found in Refs. [16-18]. Early attempts to 
calculate tri-linear coupling contributions to the neutrino masses can be found in [19-23]. A 
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basis independent calculation was performed in [24] . The complete set of 1-loop corrections 
to neutrino masses and mixings were investigated in detail in [25]. There, collections of 
weak-scale couplings that could give rise to the tri-bi maximal mixing pattern [26] were 
suggested. Such a pattern provides a successful fit to oscillation constraints on lepton 
mixing. 

The number of parameters in Eq. ( |l.5|) may be reduced by assumptions about the 
structure and origin of the supersymmetry breaking terms. We take the concrete example 
of the mSUGRA assumption 1 , where at a high energy GUT scale Mx ~ 10 16 GeV the 
scalar masses are set to be diagonal and equal to mo, all soft trilinear couplings are equal 
to their analogous superpotential parameter multiplied by a common coupling Aq and all 
gaugino masses are set to be equal to Mi/ 2 - Constraints upon such a model including LNV 
effects have been studied in Refs. [27,28]. In the latter, a cosmological bound Y2i m ^ < 0.7 
eV was placed upon the model but no attempt was made to fit the neutrino properties 
deduced from oscillation measurements. In Ref. [29], three relatively large A^- fc defined 
in a weak basis at the GUT scale were used to generate K{ and A' i3 3 at the weak scale. 
These couplings are required to satisfy upper bounds obtained from neutrino oscillation 
data. The input parameters are not directly related to neutrino phenomenology, but lead 
to interesting collider signals. However, no detailed fitting of the neutrino oscillation data 
was attempted. Bi-linear contributions to the neutrino masses in LNV mSUGRA were 
considered in Ref. [30]. Solutions with 3 non-zero GUT scale LNV parameters Hi were 
found. Also, a preference for the vacuum oscillation solution [31] to the solar neutrino 
problem was stated. This has since been ruled out by KamLAND data [32]. An attempt 
to fit the neutrino masses in LNV mSUGRA using Aj33 and A'^ in the charged lepton 
mass basis was made in Ref. [33] using the renormalization group equations (RGEs) of 
Ref. [34]. Only the atmospheric mixing angle and the neutrino mass squared differences 
were included in the fit and a preference towards the small mixing angle MSW [35] solution 
for the solar neutrinos was found. The latter has also been excluded by KamLAND, which 
favours the large angle MSW region of parameter space. 

In the present paper, we shall follow an alternative approach. At the GUT scale and 
in the weak interaction basis where all soft terms are diagonal in flavour space, we shall 
assume zero bi-linear LNV operators. In this basis, we shall introduce a small number of 
non-zero tri-linear LNV operators, following Ref. [28]. Further non-zero bi- and tri-linear 
couplings are generated on renormalisation to the weak scale. One may hope that a single 
GUT-scale tri-linear operator could then be enough to generate neutrino masses. We shall 
find that a single non-zero trilinear coupling is insufficient to reproduce the neutrino mass 
pattern observed. A further operator shall be necessary. Although we assume only two 
non-zero LNV couplings, we think of this only as a limiting case, where any other LNV 
parameters contribute negligibly to the neutrino masses and hence are neglected. Such 
flavour structure may arise from some fundamental flavour physics models which yield a 
small number of dominant LNV couplings in the weak interaction basis, mirroring the case 
of the Standard Model Yukawa couplings. 

1 Also termed the constrained MSSM. 
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The paper proceeds as follows: section [2] explains the difficulty in obtaining the ob- 
served neutrino mass hierarchy with only one R-parity violating (RPV) 2 parameter at Mx 
in any flavour basis, and suggests how the next to minimal case with two RPV parameters 
(defined in a mixed charged lepton basis) may get around the problems. We then describe 
a numerical procedure to find the best fit value for the 2 RPV parameters, as well as the 
3 mixing angles which define the mixed charged lepton basis in section ^. The results 
are presented in section f|, where we briefly speculate about possible collider signatures. 
We compare our paper with selected work in the literature, and discuss issues related to 
parameter tuning before concluding in section [|. We set our notation in Appendix [A|, and 
the general expressions for the 1-loop neutrino-neutralino mass corrections in Appendix |B]. 
Appendix [C] describes the mass insertion approximation for CP even and CP odd neutral 
scalar contributions (required for numerical stability), and finally in Appendix ^ a full 
calculation of the branching ratio of I — > I 7 is presented, which will be used to constrain 
the LNV parameters. 

2. Neutrino masses and mixings 

In a recent global 3 neutrino fit to all oscillation data [36] , the following ranges of parameters 
at la were found 

Amli = 7.9±g;|| x 10" 5 eV 2 , |Aml x | = 2.6 ± 0.2 x 10" 3 eV 2 , (2.1) 
sin 2 # 12 = 0.31 ± 0.02, sin 2 # 23 = 0.47±g;g|, sin 2 # 13 = 0^;[J 08 . 

Here Am^ and Am^ denote the mass squared differences between the three neutrinos 
responsible for solar and atmospheric oscillations repectively. The sin 2 #'s represent the 
rotation angles characterizing the PMNS matrix [37], the lepton counterpart to the CKM 
matrix, in the standard parameterization [38]. First, we shall discuss some issues concerning 
mixing matrices and effective neutrino and charged lepton mass matrices in the presence 
of lepton number violation. 

Since we work in the weak interaction basis, it is instructive to understand how diag- 
onalization of the charged lepton Yukawa matrix behaves, particularly under renormaliza- 
tion. In the R-parity conserving limit and in the absence of the neutrino Yukawa matrix, 
one can simply rotate the leptonic superfields back to a basis where all the leptonic terms 
are flavour diagonal at Mx- The rotation matrices Zn and Z[r are defined by 

Z\ L Y E Z IR = Y E , (2.2) 

where Ye is diagonal. Since mSUGRA soft terms are diagonal in flavour space, there is no 
intrinsic flavour violation in the RPC limit. It is therefore always possible to diagonalise the 
leptonic Yukawa couplings and the slepton mass matrix by the same rotations even after 
renormalization to lower scales. The rotation is also renormalization scale independent. 

To show this we examine the renormalization group (RG) equations of Ye [39] , keeping 
only terms proportional to Yukawa matrices, as the gauge coupling contributions are flavour 

2 We will use RPV and LNV interchangeably in this paper. 
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blind and so cannot change the lepton mixing. At 1-loop, we have 



dt 



Y E 



(47T) 



-Y E {T,{2,Y D Yl + Y E Yl) + 3Y E Y E + g.c.} 



= aY E \Y E Y E + $ (Tr(Y D Yl), Tt(Y e Y e ), 9l , g^j } , 



(2.3) 



where g.c. stands for gauge contributions and t = In fi is the natural logarithm of the 
DR renormalisation scale. The function p( r Tt(Y E )Y E) ),'Tr(Y E Y E ), g\, g 2 ) depends on the 



trace of the quark and lepton Yukawa matrices as well as the gauge couplings, and so is 
proportional to the unit matrix and do not change lepton flavours as a result. 
Writing Y E {t) = Y t , (3(t) = f3 t and Y E (0) = Y , we obtain 



Y t = Y + a dhY tl (Y t \Y tl +4). 
Jo 

We may then plug Eq. (|2.4j) back into itself in order to obtain an expansion 



Y t = Y + a / dt! 



J 

Jo 



Y + a I 1 dt 2 Y t2 (Y t \Y t2 + 4 
Jo 

o 



y + a/ dt 2 Y t2 (Y t \Y t2 + t2 ) 
o 



(2.4) 



(2.5) 



Eq. (2.4) may be substituted back into Eq. ( |2.5| ) an arbitrary number of times, but it is clear 
that all terms are products of terms which are either of the form f dt n Yt n (Y^Yt n )D(t n ) 
or J dt n Y tn D' (t n ) where D(t n ) and D'(t n ) are some functions of t n which are diagonal in 
flavour space. We may now calculate y to arbitrary order in a, for example 



Y t = Y + a \aY + tY Y*Y 



?Ly Y*YoYSy + BYqYJYo + CY 



+ 0(a 3 ),(2.6) 



where A = f^dhP^, B = 3 f* dtx f Q h dt 2 $t 2 + f Q t dtit 1 $ tl and C = J *cfti/3 tl f Q h dt 2 (3t 2 - It 
is clear that higher order corrections always have the form y^y^y))™ plus some diagonal 
piece, and so both Yt and Yq can be diagonalized by the same mixing matrices Zii,m\ . This 
procedure can be generalised to higher loops, but we decline to do so here. Intuitively, in a 
non-diagonal Y E basis, the mass eigenstates are defined by the vectors in the lepton mixing 
matrices Z\l and Zm, and after the rotation the new flavour basis coincides with the mass 
basis. Because there is no intrinsic lepton flavour changing physics in the Lagrangian, these 
mass eigenstates cannot mix with each other under renormalization, except through a trace 
which then affects all three flavours in the same way. The mixing matrices hence remain 
the same. 

This means that even though we work in a non-diagonal lepton flavour basis, lepton 
flavour changing processes only arise from the LNV couplings and are suppressed by even 
powers of them since the LNV operators violate lepton number by one unit and the Ma- 
jorana neutrino masses generated violate it by two. Another source of leptonic flavour 
changing neutral current (FCNC) would arise when one allows for non-universal SUSY 
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breaking boundary conditions, which is common in unified model building. This is a sep- 
arate issue and has been discussed extensively in the literature, see e.g. [40]. 

The LNV operators also make the rotation matrices Zn and Zm scale dependent. 
However as R-parity violation is small in our analysis, the mixing matrices are scale inde- 
pendent to a very good approximation. For numerical convenience, we rotate to a diagonal 
charged lepton Yukawa matrix basis at Mx before renormalizing down to Mz- In general 
this means that there are a lot of LNV parameters in the diagonal basis, however they are 
all related to each other by the rotation. 

We will first investigate the case with one non-zero dimensionless LNV coupling at 
Mx, and then go on to discuss the case of 2 non-zero LNV parameters. The RG evolution 
will dynamically generate all LNV parameters allowed by symmetry. This includes the 
dimensionful /Vs, which implies non-zero sneutrino vacuum expectation values (vev) upon 
minimizing the Higgs-sneutrino potential [28]. These effects are important in determining 
the neutrino mass pattern and will be described in detail in the following sections. The 
mass matrices used to compute the radiative corrections of the neutrinos and neutralinos 
are presented in Appendix [A] and references therein. 

The tree level, 7x7 neutrino-neutralino mass matrix M. x generates an effective 3x3 
neutrino mass matrix through a seesaw mechanism. Writing A4 x as defined in Appendix [A|, 

M N = * , (2.7) 

V m m v J 

where m u and m contain the lepton number violating contributions, and M x o is the neu- 
tralino mass matrix which plays the suppression role of the right handed heavy neutrinos 
in the standard see-saw mechanism [7]. The effective mass matrix A4g ff is then given by 

A< ff = m u - mM7,lm T . (2.8) 

At tree level, m v = 0^x3 and the effective mass matrix is given by [9, 10] 

/ ApAp A fi A„ A«A 



M 



V 

off 



(M l52 2 + M 2 g 2 



2n[v u v d (M ig % + M 2 g 2 ) - l iM 1 M 2 ] 



A M A e A^A M A^A r | . (2.9) 
y A T A e A r A M A T A T 



g 2 ,g are the SU(2) and U(l)y gauge couplings, and v u ,Vd are the vacuum expectation 
values (vevs) of the neutral components of H u ,Hd respectively. Also, 

Ai = nvi - v d fM, i = {e,fj,,r}, (2.10) 

and so depends on the sneutrino vevs Vi and LNV supersymmetric bilinear parameters 
fj,i. is of rank 1, so there is one non-zero eigenvalue 

(M lg j + M 2ff 2)EA 2 
mheavy 2n[v u v d (M l9 l + M 2 g*) - vM x M 2 \ ' ^ 

Clearly this is not compatible with the experimentally measured mass squares differ- 
ences, which require at least two non-zero neutrino masses. Including one-loop contribu- 
tions can give mass to the other neutrinos. In the rest of this work we will present a 
minimal construction in which the observed mass pattern can be obtained. 
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The complete one-loop self-energies for the neutralinos and neutrinos in 2 component 
notation, including all terms with R-parity violation are given in [25]. Here we present the 
general expression in a matrix form compatible with the notation of [41]. The one-loop 
corrected neutrino-neutralino mass matrix M' N is given by 

(M' N )ij = {M N )ij + \{&M + 6M T )ij, (2.12) 

where 



SMij = & D )ij ~ {M N ) ik (X L ) kj , (2.13) 

and £d and T,l are mass corrections and wavefunction renormalization respectively. 

As in [41], we perform the full renormalization with external legs in a flavour basis. 
The Feynman rules including all LNV terms can be obtained in the Appendix of [25], 
and by rotating one external neutral fermion leg from mass basis into flavour basis. The 
physical masses are obtained by diagonalising M.' N . 

Upon diagonalising M' N , a 7 x 7 mixing matrix Z' N is obtained. As at tree level, one 
may also obtain an effective neutrino mass matrix -M^ using Eq. ( [2.8D . In particular, 
the 3x3 sub-matrix m u receives radiative contributions and plays a significant role in 
determining the effective neutrino mass matrix. The mixing between the neutrinos and 
the neutralinos is suppressed by the smallness of the LNV couplings and Vi with respect 
to the other neutralino matrix parameters, and is neglected for the purpose of calculating 
the PMNS matrix. The neutrino mixing matrix Z v is then defined by 

Z^M%Z U = diag[mi], i = 1, 2, 3. (2.14) 

The charged lepton-chargino mass matrix is treated similarly, where one can again use 
a see-saw mechanism to obtain an effective mass matrix with the charginos providing the 
see-saw suppression. In this case however, the see saw contribution to the effective mass 
matrix is suppressed by two powers of LNV parameters, and is negligible compared with 
Ye- So for the purpose of calculating the charged lepton mixing matrix, one can simply 
use Ye and hence the same Zu, and Zm as before. After rotating to a basis with diagonal 
lepton Yukawa couplings at Mj, the lepton flavour changing terms in the Yukawa matrix 
Ye are generated by terms at least quadratic in A by the RG evolution, and again such 
contributions are neglected when calculating the PMNS matrix. 

The observable PMNS mixing matrix which connects the charged lepton and neutrino 
mass eigenstates is defined by 3 

Upmns = Zf L Z v , (2-15) 



3 Note in 4 component Dirac notation, the Ye is typically associated with the term LLYsEphd in other 
conventions. This Ye is the complex conjugate of the Ye that appears in the superpotential in Eq. ( ]l.3| ) . This 
is why we have Upmns = Zj L Z v instead of Upmns = Z\ L Z V . However our numerical calculation assumes 
CP conservation, so the two are equivalent. The results presented in this paper follow the convention 



defined by Eq. (1.3), and we keep track of the complex conjugations for book keeping purposes. 
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where Zn, is the mixing matrix of the left handed charged leptons. In the present paper, 
we always work in the CP-conserving limit and thus z\ L = Z^ L . Because we are already 
working in a basis where the charged lepton mass matrix is practically diagonal, the PMNS 
matrix is the same as Z v , and they will be used inter-changeably in the rest of the paper. 

The ratio of tree level to loop level mass scales 

It is well known that even if the high scale theory has all bilinear LNV couplings set to 
zero, the RG flow will dynamically generate /ij, Di and m\ dL . at low scale, and hence «j 
after minimizing the higgs-sneutrino scalar potential. These effects are important, because 
the dimensionless LNV couplings contribute to the effective neutrino mass matrix at both 
loop level and tree level in a related way. Typically the tree level mass scale dominates 
over the loop contributions, despite the fact that the bilinear parameters themselves are 
originated from radiative corrections. Also, the loop induced mass matrix tends to align 
with (i.e. is proportional to) the tree level mass matrix, further suppressing radiatively 
induced neutrino masses. 

With only one trilinear LNV pa- 
rameter in the diagonal charged lepton 
basis, the above statements are gener- 
ally true. However one LNV parame- 
ter in the weak interaction basis cor- 
responds to many LNV couplings in a 
diagonal charged lepton basis, related Figure 1: Dynamical generation of ^. 

by the charged lepton mixings. These 

rotations need not be small because the observed PMNS mixings are near tri-bi maximal, so 
we may expect a 'democratic' distribution of LNV couplings. However the fact that these 
parameters are related implies that both the alignment effect and the hierarchy between 
the tree and loop mass scale persists. We will show how this happens in the following, and 
then go on to discuss the next minimal case of two LNV operators (in the weak interaction 
basis) . 

In scenarios with R-parity violation dominated by one A, one can approximate the tree 
level mass scale by [28] (see also [42] ) 4 

m% ^ -!^[^ ?]2[ i n ^ ]2[ A ii9( F|) i? ]V 2 (4;4;4;t-/3), (2.16) 

mi mi mi 




5M 1/2 T6vr 2J L M z 



1/2 -L07T- iv±z mjj mj, 

where / is a dimensionless function of O(10). To understand the form of this expression, 
note that the tree level neutrino mass is obtained by Aj in Eq. ( p.lOj ), which is proportional 
to Hi and the sneutrino vevs Uj, while the latter also scale like /i, up to effects from 
the Yukawa couplings. The dynamical generation of /Xj, neglecting sub-dominant terms 
proportional to ^ itself, is determined by the RG flow 

16tt 2 ^ = -n(*i jk (Y£) jk + 3% jk (YZ) jk ), (2.17) 



'Note that the normalization of the higgs vev in the current paper is different from [28]. See Appendix^. 
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so a first estimate of \ii at Mz may be obtained by naively integrating the right hand side 
of the above equation, yielding 



* = (4^M^)(A«*(yS)i* + 3^*05)^ + ((^4^^) • (2-18) 



A pseudo-Feynman diagram in terms of superfields for the second term of Eq. ( 2.18j ) is 
shown in Fig. [l|. 

On the other hand, the scale of the radiative corrections can be approximated by the 
part of T,£) in Eq. ( f2.13| ) which correspond to m v in Eq. ( |2 .8| ) (see also [25]). We define 
this scale to be 

^ m * {M f ) kk (M) ) kk 

loop = \^ /y. \ \ ..,\., ll StR It, L" 



'"" s £ M - E A w% fg (2 ' 19) 

i=l ij.fc III }rr Irr 



where the summation i is over external neutrino flavour eigenstates in the basis defined 
in Appendix [A]. m*j 5 is the chirality flip of the fermion propagator with legs labelled by 
j, and the M 2, s represent mass insertions with appropriate scalar 'handedness'. A mass 
insertion diagram of m° op is displayed in Fig. ^. The approximation for A' is exactly the 
same, apart from extra colour factors n c = 3 for both Eq. ( 2.16| ) and Eq. ( |2.19|) . 
The ratio mfj' ee /mi oop can then be approximated to be 



mi ree 



a GUT m\M x /M z ) (M}J - QMy ^ ^ j_ 
loop- c B7rM n ,nMn - ;/.ta.Ti/?l Inf7/W 2 - n v 



5vrM 1/2 (A)-Mtan/3)ln((>K )/(M 2 ? )) J X m 2 ' m 2 

JLL J RR 



where we have made the approximation 



(M 2 hR ) kk ~ m /fc (^ - Man/3), (2.21) 



where / represents sleptons or squarks, as appropriate. Setting aci/T = 0.041 and Mx 

h/2 



10 16 GeV, a rough scan along the SPS1 slope [43] with M 1/2 = [250,750] and also tan/? 



[10,40] gives mi ree /ml oop of order 0(30 - 80) for A^-, and 0(70-200) for A^-, in agreement 



with the approximation of Eq. ( 2.20| ). However, this is too large to account for the mild 



empirical neutrino mass hierarchy of a factor of <6 from Eq. ([2.1|), assuming negligible 
lightest neutrino mass 6 . 

Note that in the above discussion, alignment (proportionality) between the tree level 
neutrino effective mass matrix and the radiative correction has not yet been included. 
The inclusion of such effects will enhance the hierarchy between the heavy and the light 
neutrino masses obtained in the above approximation. Because at tree level there is only 
one massive neutrino, exact alignment would imply still only one massive neutrino after 
the inclusion of radiative corrections. Misalignment between the tree level contributions 



5 The phases of the charged lepton mass matrix rrif can always be rotated away by redefinition of 
the lepton fields. However we find leaving the redundant complex conjugations helpful when considering 
chirality flips in loop diagrams. 

6 In the opposite extreme of a quasi-degenerate mass spectrum, there is clearly no mass 'hierarchy'. 
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and the radiative corrections will give rise to light neutrino masses. The size of the mass 
prediction for the second lightest neutrino compared to the heaviest neutrino, and hence 
|Am|i| compared to |Am§i|, can be estimated by the amount of misalignment, multiplied 
by the ratio of the loop/tree level scales. 

To discuss possible alignment effects, / . >o* \ 

it is instructive to first start in the charged 

lepton mass basis with one LNV cou- Vi ► 

pling. In the latter case, if the LNV cou- *ijk V J ^ikj 

pling violates only one lepton flavour, m* f 
then the other two lepton flavours are 

separately conserved. Each of the tree Figure 2: A mass insertion diagram which repre- 
level and loop level neutrino mass ma- sents the loop corrections to m l ° op by Ay*, and \kj ■ 
trices can only have one non-zero entry 

corresponding to the lepton flavour violated, and so these contributions to the full neutrino 
mass matrix are exactly aligned. Using this flavour symmetry argument it is easy to see 
that this alignment will persist to all orders in perturbation theory. On the other hand, 
if the LNV coupling violates all three lepton flavours, then no bilinear LNV couplings can 
be dynamically generated. There is no source which violates lepton flavour by two units, 
therefore we cannot generate a one-flavour violating bilinear term from a three-flavour vi- 
olating term. This can be seen explicitly in the one-loop approximation by reference to 
Eq. ( p. 17 ). Thus no tree level mass matrix can be generated. The same reasoning, together 
with the requirement of an odd number of chirality flips in the loop diagrams, also implies 
that the radiative correction to the neutrino mass matrix is zero. 

When the one LNV coupling is defined in the weak interaction basis, which need not 
coincide with the diagonal charged lepton basis, exact alignment will cease to hold, as the 
above symmetry arguments do not apply any more. However, it is possible to show that 
to lowest order in perturbation theory, there is still alignment between the tree level and 
loop level contributions. For this purpose, we work in a diagonal lepton basis, and make 
use of the fact that in the R-parity conserving limit, the mixing matrices that rotate into 
the diagonal basis obey 

Z lL = Z lL Z lR = Z~ m (MlJ jk cx(Y E ) jk . (2.22) 

Eq. ( 2.22] ) is true at all renormalization scales, as discussed in the previous section with 
universal boundary conditions due to the absence of intrinsic lepton flavour violation. The 
universality of the lepton and slepton mixings will be violated at order A 2 , but we shall 
study small values (< 10 -3 ) of A and so the non-universality is neglected in our lowest 
order approximation. We will also make use of the fact that all dynamically generated 
LNV couplings are small, as will be seen shortly. 

In the following we denote parameters in the diagonal lepton basis with a hat. Although 
we use A for illustration, the case of A' is similar. The dimensionless LNV couplings in the 
diagonal lepton basis, \ a bc are related to the LNV couplings in the weak interaction basis, 
^ijk, by 

^abc = \jk{Z\ L )ai{Z\ L )bj{Zm)kc- (2.23) 
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The tree level neutrino mass matrix then scales (neglecting the sneutrino vevs) as 



and so 



{■Mtfdab oc A a A b oc fi a fi b oc A acd (y|) cd A te/ (y|) e/ 

nb ■ 



(2.24) 



(2.25) 



where M^ s is the effective neutrino mass matrix in the original lepton flavour basis. The 
factorization of the two charged lepton mixings represents the shift of Zil from the charged 
lepton sector to the neutrino sector when computing Upmns using Eq. ( p.!5| ). At loop 
level, we approximate its contribution to the full mass matrix by the dominating effect 
from E^> 



ab 



EE 

c,f d,e 



[m 



l)cf 



(4vr; 



2 X acd \ be fB o (0,m~ ■ ,™>i af ) + (a^b) 



(2.26) 



in self-evident notation, and where Bq is a Passarino-Veltman function [44]. Using a mass 
insertion approximation (MIA), one obtains 
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+(t «-» n) }(z} L ) ai (Z? L ) nb , 



(Mf ) dd 

l RR 



)] 



1 + 



(2.27) 



where -Mr is a mean squared mass of the left (right) handed sleptons, and 5M 2 rep- 
resents deviation from such mean values. Note the index contraction of the A in Eq. ( 2.27| ) 
is different from those in Eq. fl2.24| ). However is dominated by the one LNV coupling 
that is non-zero at Mx, as all dynamically generated couplings are small compared with 
it. We then have i = n, j = I and k = m and consequently Ti v D is proportional to the tree 
level effective mass matrix M^ s , i.e. 



(±» D ) ab oc (M^) ab , 



(2.28) 



up to higher order corrections. It is unwieldy to write down the higher order corrections to 
this expression, with terms from many sources (e.g. RGE, mixing matrices including the 
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charginos-lepton mixings and corrections in MIA) entering the next lowest order pertur- 
bation. The key result here is that the index contractions in Eq. (2.24) and Eq. ( 2.27| ) are 



different, unless they are dominated by a single LNV parameter, in which case they become 
proportional to each other. The tree level mass matrix is bilinear in Aj, which is propor- 
tional to /ij. The dynamical generation of \i{ can be obtained by integrating their RGEs, 
which contain terms with odd powers of Xijk and A^- fc . On the other hand, one obtains the 
loop corrections to the neutrino mass matrix by summing over contributions which contain 
even powers of the LNV couplings. The interference between various contributions at tree 
level and at loop level are in general different as a result. 

The scale of the light neutrino mass can be estimated by noting that in the case of 
degenerate charged lepton masses, there is an ambiguity in the definition of the charged 
lepton mixing. The Yukawa matrix is simply a unit matrix, which implies that our earlier 
discussion on one LNV coupling in a diagonal charged lepton basis will apply in this 
situation, and so there is only one massive neutrino. In other words, the light neutrino 
mass m]] ght arises from the deviation of charged lepton universality. Its mass is therefore 
expected to be suppressed relative to the heavy neutrino mass m^ eavy by (the square of) 
the ratio of the charged lepton mass to the scalar mass scale that appears in the loop, as 
well as the tree-loop scale difference coming from the RG evolution, ie 

- 0(&'\ (2.29) 



heavy \ m~ 



with / = A; tan/3) the function from Eq. (|2.20|) 



p "'0 

With more than one dominant LNV coupling, the situation is very different. By 
choosing appropriate LNV parameters, one can obtain a partial cancellation between the 
dimensionless LNV contributions in the dynamical generation of /ij, and suppress the tree 
level neutrino mass as a result. Also the above derivation of alignment relies on one 
coupling dominance, specifically i = n, j = I and k = m Eq. ( 2.2^ ). Once we have more 



than one dominant coupling, this manipulation ceases to hold. As the alignment effect is 
weakened and the mass hierarchy becomes smaller, larger PMNS mixing is expected with 
large charged lepton mixing angles. So in the next minimal case with 2 LNV couplings 
at Mx, by allowing partial cancellations in the fii generation and large charged lepton 
mixings, we may already be able to obtain the observed neutrino mass pattern. We will 
investigate this possibility by a numerical approach, discussed in the next section. 

3. Numerical procedure 

The parameter set that defines our model is given by 

mo,M 1/2 , A),sgn//, @ M x , 

e[ 2 ,e[ 3 ,e l 23 ,A 1 ,A 2 @ m x , 

tan/3 @ M z , (3.1) 

where 

Ai,A 2 G {\i jk , ( 3 - 2 ) 
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and Mx is denned to be the scale at which the electroweak gauge couplings unify, mo, 
M1/2 an d Ao are the universal scalar mass, gaugino mass and SUSY breaking trilinear 
scalar coupling (divided by the corresponding Yukawa coupling) at Mx respectively. The 
sign of [i is defined by sgn^x. The charged lepton rotations are characterized by 6 l, s in 
the standard parameterization [38], and tan/3 is the ratio of the higgs vacuum expectation 
values v u /vd- 

In the quark sector, flavour mixing is assumed to reside in the down-quark sector at the 
weak scale. We have also assumed that the down-quark Yukawa matrix Yd is symmetric 
at that scale. In the lepton sector, we work with non-diagonal charged lepton Yukawa 
matrices, and define any non-zero LNV parameters in that flavour basis. The mixing 
angles 9 l are determined by fitting the PMNS matrix. For definiteness, we assume that the 
two charged lepton mixings Zn and Zm are the same, or equivalently, that the charged 
lepton Yukawa matrix Ye is symmetric. 

The set of Lagrangian parameters (RPC and LNV) at Mz are obtained using a modified 
version of the SOFTSUSY code [45]. At Mx, a set of dimensionless LNV parameters defined 
in Eq. ( |3.1| ) is chosen. At Mz, tan/3, the ratio of higgs vevs is also specified. The high scale 
LNV parameters are rotated to a diagonal charged lepton basis specified by 9 . 

The parameters are run to Mz using the full set of 1-loop MSSM RGEs including 
all LNV contributions, which are then required to match with the low scale boundary 
conditions on the gauge and Yukawa couplings. Throughout the analysis, we use nit = 
172.7 GeV [46] and Mz = 91.1876 GeV for the pole masses of the top quark and the 
Z° boson respectively. The weak scale gauge couplings in the MS scheme are set to be 
a~ l {Mz) = 127.918 and a s {Mz) = 0.1187. Light quark masses are also set to their central 
values in the MS scheme: mb(mb) = 4.25 GeV, m c (m c ) = 1.2 GeV, m s (2GeV) = 0.1175 
GeV, m d (2GeV) = 0.00675 GeV and m u (2GeV) = 0.003 GeV [38]. For the calculation of 
the physical Higgs mass and electroweak symmetry breaking (EWSB) conditions, one loop 
corrections are included. 

To calculate the neutrino masses, we run the parameter set to 

M S = yJm h (M s )m h (M s ), (3.3) 

where the EWSB conditions are imposed and all sparticle and Higgs pole masses are calcu- 
lated in SOFTSUSY. The neutrino masses are obtained from the full 7x7 neutrino-neutralino 
mass matrix including all LNV contributions at the 1-loop level. Throughout the calcula- 
tion, the (s)lepton doublets are put on the same footing as Hj. Our analytic results were 
checked against Ref. [25] for the full LNV case, and against Ref. [41] in the RPC limit. Our 
numerical results agree well with Ref. [25]. In the RPC limit, the neutralino corrections 
are also in agreement with the latest published version of SOFTSUSY. 

After the initial rotation to a diagonal Ye at Mx, the charged lepton matrix remains 
essentially diagonal at all renormalization scales due to the smalless of the LNV couplings. 
The computation of the PMNS mixings thus neglects all 'residual' flavour mixing generated 
through renormalization in the charged lepton sector. The three charged lepton mixing 
angles, together with the two GUT-scale LNV parameters are varied in order to obtain the 
best fit. 
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The minimization to find the best fit parameter set is performed using MINUIT [47] 
and the observables assumed in Eq. ( |2.1| ). x 2 1S defined to be 

-^obs / £ s \ ' 

2 _ / Ji ~ 



where Oi are the central values of the N Q b s experimental observables, fi are the correspond- 
ing numerical predictions, and <7j are the 1-sigma uncertainties. We find that it is difficult 
to obtain convergence of a procedure which minimizes x 2 when solar and atmospheric mass 
squared differences are included in Oi directly. We find a multi-step minimisation approach 
to be effective. In this approach, we first minimize a x 2 °f t ne mass ratio of the 2 heavy 
neutrinos and the sine squared of the three PMNS mixing angles by keeping one LNV 
parameter fixed. For a normal hierarchy we have 

m Ul ~ 0, m U2 = 8.9±{j;i6 x 10~ 3 eV, m Vz = 5.10 ± 0.20 x 10- 2 eV, (3.5) 

ma 

=> — = 5.74 ±0.32, 
m 2 

whereas for an inverted hierarchy, we have 

m Ul = 5.10 ±0.20 x 10~ 2 eV, m V2 = 5.18 ± 0.19 x 10~ 2 eV, m U3 ~ 0, (3.6) 

^ = 1.0151 ±0.0769. 

mi 

At this stage, we obtain a good fit to the charged lepton mixing angles and the mass 
ratio but the values of the neutrino masses in general do not fit the atmospheric and solar 
mass squared differences. To perform a 'proper' fit we then switch to a x 2 °f Am 2 ,-^ l^ m ii I 
as well as the three PMNS angles. We first fix the three charged lepton mixing angles and 
allow the two LNV parameters to vary. The LNV parameters will then adjust their values 
to fit the mass squared differences. Finally a full five parameter fit is performed. It is also 
helpful to have a good initial estimate of the size of the LNV parameter that is fixed in the 
first stage of minimization. We find that a suitable starting point is the value of the LNV 
coupling that saturates the one LNV coupling bound in [28] , derived from the cosmological 
bound Y^i m u z < 0.7 eV at 95% condidence level (CL) [48]. 

The above fitting is performed at Ms- After this, the resulting set of parameters is 
renormalized down to Mz, where specific products of couplings are checked against the 
bounds in [27]. We also checked against the experimental bounds on the branching ratios 
of the processes I — > I'j [49-51], B s — > [52] and b — > S7 [53]. The bounds we use are 



BRexpiH ■ 


-> erf) 


< 1.2 x 10- 11 


90%CL 


BR exp (T - 




< 1.1 x 10~ 7 


90%CL 


BR exp (r - 


-* m*) 


< 6.8 x 10~ 8 


90%CL 


BR eX p{B s - 




) < 1.0 x 10~ 7 


95%CL 


BR(b -i 


■ «7) 


< 4.34 x 10~ 4 


2a, 



(3.7) 

BR exp (B s -> n + n~) < 1.0 x 10~ 7 95%CL, 

2.76 x 10~ 4 < 
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where BR(b — * 57) combines statistical and systematic experimental errors [54] together 
with a Standard Model theoretical error [55] by adding them in quadrature. 

We have included a detailed calculation of BR(l — > I'j) in Appendix[TJ. For a discussion 
of contributions of different couplings to this process, see e.g. [56]. The LNV contribution 
to BR{B S — > /x + /i~) was calculated in [57]. A generalisation to include left-right mixings 
of the mediating squarks and non-degenerate sparticle masses can be found in [53]. The 
LNV contribution to BR(b — > 57) in a leading log approximation is presented in [58]. 

One technical aspect that should be mentioned concerning the numerical calculation is 
the cancellation between the CP even (CPE) and the CP odd (CPO) scalar contributions 
to the neutrino masses at loop level. In the R-parity conserving limit, the CPE and CPO 
sneutrino contributions cancel exactly. When R-parity is violated, the higgses mix with the 
sneutrinos through the sneutrino vevs and the bilinear soft terms Di, lifting the degeneracy 
of the masses of the CPE and CPO scalars. The mixing matrices of the CPE and CPO 
scalars also differ. This degeneracy lifting is in general very small compared with the 
sneutrino mass scale, but may have a significant impact on the neutrino masses. 

Such a large cancellation is numerically unstable. In the phenomenologically interesting 
LNV parameter space, numerical fluctuations caused by such an instability are actually 
insignificant for weak scale sneutrinos. However away from this 'good' LNV region, the 
second-lightest neutrino mass is usually highly suppressed compared to the third lightest (or 
the heaviest) neutrino as discussed in the previous section. Its eigenvalue can be dominated 
by numerical fluctuations, and tends to spoil the minimization procedure when the LNV 
parameter space away from the physically interesting region is accessed. 

To get around this problem, we perform an analytic approximation commonly known 
as mass insertion approximation (MIA). We apply this technique to calculate the devia- 
tion from exact cancellation in the RPC limit, instead of performing the large numerical 
cancellation directly. In this calculation, it is useful to think of the 5x5 neutral scalar 
matrices as a 2 x 2 higgs block, a 3 x 3 sneutrino block, and a 2 x 3 block that mixes the 
higgses and the sneutrinos in the basis (h u , hd, v e , u^, D T ) as follows 

f M 2 \ _ (( M2 H(A)hx2 (^(A))2x 3 \ 

[M CPE{CPO) )^ - y K(a))sx2 (M|)3x3 ) • (3-8) 

We take advantage of the fact that in the 't Hooft-Feynman gauge, the CPE and CPO 
sneutrino blocks M~ are identical and can be diagonalised by the same orthogonal ma- 
trix. The difference between the CPE and CPO contributions then come from Mjj^ and 
<J H(A)- The CPE and CPO M 2 's can be separately diagonalised, leaving the off-diagonal 
contributions entirely in the (s)lepton number violating sector. In this basis it is possible 
to obtain successive corrections to the mixing matrices as well as the mass eigenvalues, 
or in our case their difference in terms of the small LNV parameters. We refer interested 
readers to [59] for a review of MIA. For details of our specific calculation, see Appendix |C]. 

4. Results 

For concreteness, we confine ourselves to the SPSla [43] point with non-zero LNV couplings, 
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i.e. 



M 1/2 = 250GeV, 
sgn/x = +1, 



m = lOOGeV, 
tan/3 = 10, 



An 



-lOOGeV, 



(4.1) 



although the numerical procedure can be applied equally well in other mSUGRA region. 
We first display an example with a normal hierarchy with non-zero A233 and A 233 at the 
The set of best fit parameters is given by 9 l l2 = 0.460, l 13 = 0.389, 



unification scale. 
0.305, A 233 



7 23 



4.07 x 10" 5 and A' 



233 



-2.50 x 10 . The variation of the observables 



charaterising the neutrino mass pattern with A 233 is displayed in Fig. ||. We see in Fig. 3(a) 



that as A233 increases, the mass difference between the two heavy neutrinos decreases due 
to suppression of the tree level mass matrix and a weakening of the alignment effect. 
Eventually the mass difference increases again, signalling a cross over where the heaviest 
neutrino becomes dominated by contributions derived from A 233 . In between there is a 
window which allows the mass squared differences to fall within the experimentally observed 
range, as displayed in Fig. 3(b). The variations of sin 2 #i 2 and sin 2 #i 3 shown in Fig. 3[c 



and Fig. 3(d) are associated with the cross over described above, when the two mass 
eigenvectors determined by the two LNV parameters switch. 
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Figure 3: Variation of the observables characterising the neutrino mass pattern with A233 at 
SPSla in a normal hierarchy. The best fit parameters are 6[ 2 = 0.460, 9 l 13 — 0.389, 6*33 = 0.305, 
A233 = 4.07 x 10 -5 and A 233 = —2.50 x 10 ~ 6 . The strips display the empirical la bands. 

In Fig. I we show an example with an inverted hierarchy, again with non-zero GUT- 
scale A233 and A 233 . The mass squared differences is much more sensitive to the variation 
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of A233, because of the larger cancellations required to obtain two near degenerate massive 
neutrinos. This reflects that fine tuning between the a priori unrelated A233 and A 233 
couplings is required to generate the two quasi-degenerate mass scales in the inverted 
hierarchy. On the other hand, we see in Fig. 4(c) and Fig. 4(d) that sin 2 # 13 and sin 2 # 2 3 



are insensitive to changes in A 233 . This is because the lightest neutrino is almost massless, 
hence its flavour composition depends only weakly on variations in A 233 once the mass 
degeneracy of the two initially light neutrinos are resolved by the presence of A 233 . With 
an inverted hierarchy, this implies that the third column of Z v is insensitive to the changes 
in the LNV couplings, which in the standard parameterization is determined fully by #13 
and 023 7 • Of course, sin # 12 will still change as shown in Fig. [4(d)| , which again is associated 
with the cross-over of the two heavy mass eigenstates. 

For an inverted hierarchy, in principle there may be sizeable two loop corrections 
which could affect the delicate balance among terms in the effective neutrino mass matrix. 
However, the two loop corrections in question are not in the literature. We expect the 
main corrections to come from the strong interaction, and possibly from heavy flavours. 
Smaller corrections may be found, for example, in [60]. We attempted to estimate these 
2-loop effects by scaling the 1-loop quark contributions. For the parameter sets displayed 
in Tables [l] and ^, it was shown that an universal 10% variation leads to a change in x 2 
of the order 30 for an inverted hierarchy. However, after re-fitting all x 2 ' s g° back to the 
level before scaling, and the change in the best fit parameters is negligible. For a normal 
hierarchy, the change in \ 2 is of the order 0.01, and the changes in re-fitted parameters 
as well as x 2 are again negligibly small. Therefore, although more careful analysis of the 
two loop effects is probably required, we do not expect the possibility of finding modified 
best fit solutions with two nearly degenerate mass eigenstates in an inverted hierarchy to 
be affected. 

We have described in the previous sections that the tree level mass needs to be sup- 
pressed relative to a 'natural' value of order given by Eq. ( 2.16j ). As shown in Fig. ^ and 



Fig. ||, this can be achieved by a partial cancellation between the two LNV couplings. To 
quantify the amount of tuning required, we follow Ref. [61] and define a fine tuning measure 
A FT by 

AFT ~ 1 01n A' (42) 



where O is either Am|i or |Am 31 |, whichever gives a larger Apr- In can be seen in Fig. 3(b) 
that for a normal hierarchy, the two observables result in similar A^t- However for the 
inverted hierarchy, Apr is determined by Am^, see Fig. 4(b) . The second case can be 



understood from the fact that Am 21 arises from the difference between mi and m 2 , both 
of which are of order 0(^/|Am 31 |). The variation of Am 2l is then expected to be of order 
0(| Am^ |), the latter of which is significantly larger than the solar mass difference. We 



7 There is also a dependence on the Dirac CP violating phase 8, but it is suppressed by the smallness of 
sin 2 #13, and we assume CP conservation throughout our numerical analysis. 
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Figure 4: Variation of the observables characterising the neutrino mass pattern with A 233 at 
SPSla with an inverted hierarchy. The best fit parameters are 0\ 2 — f. 55779, 6[ 3 — 0.815115, 
9 l 23 = 0.146126, A 233 = 1.36 23 x 10~ 4 and A 233 = -5.69116 x 10 ~ 6 . The strips display the 
empirical la bands. The lines 9 13 and 9 23 are essentially flat. 



obtain numerical values for the fine tuning measure by 

MO(\' a ))-HO(x' b )) 

T ^ 1 l^hrA] (43) 

where X' b is the best fit value, and is given by 

A; = Aft x (1 - 2.0 x 1(T 4 ). (4.4) 

For the normal hierarchy, A ft works out to be of order O(10), so the level of fine tuning is 
mild compared with other SUSY parameters (see e.g. Ref. [62]). However for the inverted 
hierarchy, A ft is typically of order O(700 — 1000). In this sense, the present study shows 
a preference towards a normal hierarchy. 

In Tables |l] and we present a collection of parameter sets that minimizes the x 2 
defined in section |3[ We also display the values of Aft- Naturally not all parameter 
combinations result in good fits of the observed parameters listed in Eq. ( p.l[) . Note however 
that we do not attempt to exhaust all possibilities within this 2 LNV parameters setting. 
Also, not all LNV parameter combinations result in satisfactory numerical convergence. 
While some is due to flavour structure incompatible with the requirement of tree level 
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suppression, others are due to problems encountered in the minimization, for example 
running into a corner in the parameter space. Bearing this in mind, a few comments can 
be made: 

• Because the LNV couplings involved are all very small, the sparticle mass spectrum 
remains essentially the same as in SPSla, and all branching ratios of the sparticles 
are unaffected except for that of the LSP. 

• It is not possible to obtain the observed mass hierarchy with 2 A"s, because in order 
to suppress the tree-loop mass scale difference, the two A"s need to violate the same 
lepton flavour in order to contribute to the same in for the cancellation to occur. 
However, with only one flavour violated, there can only be one massive neutrino. 

• The successful parameter sets that fit the data all consist of one A and one A' coupling. 
In principle it is possible to have two A's, which might be more interesting because 
this implies the LSP will decay almost exclusively via leptonic channels. Indeed if 
the assumption of Zil = Zm is relaxed, we find parameter sets with two dominant A 
couplings. 

• Typically all three charged lepton mixing angles are large. This means that we expect 
many competing decay channels for the LSP, even though we only have 2 non-zero 
LNV couplings in the weak interaction basis. As an example, we display in Tables || 
and |] the major LSP decay channels with A233 and A233, for a normal hierarchy and 
for an inverted hierarchy respectively. 

• Although we focus on the SPSla point, in which the LSP is a neutralino, it is also 
interesting to consider scenarios with e.g. stau LSP. In the latter case, a back of 
an envelope calculation shows that for a 4 body decay to dominate, one typically 
requires (in the diagonal charged lepton basis) LNV couplings with no stau index to 
be roughly 4 orders of magnitude larger than the couplings with a stau index. This 
seems unlikely given the necessarily large charged lepton mixings. 

• The lightest neutrino mass is negligible in the current set up. Given this, the other two 
neutrino masses can be obtained using the atmospheric and solar mass squared values. 
Regardless of whether one assumes the inverted or natural hierarchy, the resulting 
neutrino mass sum automatically satisfies the cosmological bound Yli "V, < 0.7 eV. 
However, the size of the couplings at Mz can still be comparable to the single coupling 
bounds obtained in [28] based on the cosmological bound. This is because the required 
cancellation for the tree level mass suppression makes the couplings larger than what 
would be naively expected. 

• Roughly speaking, for sparticle mass of 0(100) GeV, the LNV coupling products 
involving AA or A'A' need to be of order ~ 0(1O -6 ) to have any chance to saturate 
the BR(fj, — > e'j) bound, and are less stringent for the other two leptonic FCNC. 
These can be estimated readily using the analytical expressions in Appendix For 
BR(B S — > fi + [i~), the coupling products need to be at least of O(10 -5 ) [57], whereas 
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for BR(b — > S7) the coupling products need to be of 0(1CT 2 ) [58]. Clearly the 
couplings involved in our cases are too small to contribute significantly to these rare 
decays, but we have also checked numerically that this is indeed the case. 

If the LSP is a stau, it is expected to decay promptly at the interaction point. If the 
LSP is a neutralino, it decays into 3 SM fermions. Assuming a pure photino LSP, an 
order of magnitude estimate for the 3 body decay is given by [63] 



r(x? - fff) 



n c a(X or A 
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(4.5) 



where the /'s denote SM fermions, n c is the number of colours, a is the fine structure 
constant, and Msusy is the mass scale of the virtual, intermediate scalar. For the 
results displayed in Tables [l] and E| the neutralino can decay with a displaced vertex 
of O(0.1)mm (the total width is of order 0(1O -12 — 10 _11 )GeV). Such a displacement 
would not be immediately obvious, but could be searched for, providing additional 
confirmation. 

For the inverted hierarchy, due to the high degree of tuning, all digits presented in 
Table are needed in order to produce the correct neutrino oscillation parameters. 
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= 4.09870 x 10" 


5 


2.01028 


1.07090 


2.18984 


8.00 


%n = 


-1.14183 x 10" 


-3 


^321 


= 9.88004 x 10" 


5 


0.448302 


0.400044 


1.32004 


12.4 


f A 3 ii = 


-8.66312 x 10" 


-4 


Al32 


= 9.85792 x 10" 


5 


0.921694 


-1.12754 


2.00207 


9.26 


'A233 — 


-4.07409 x 10" 


-5 


A 211 


= 8.36075 x 10" 


4 


0.457859 


0.390978 


0.304896 


8.08 


^233 — 


-4.07226 x 10" 


-5 


\l 

A 221 


= 4.16596 x 10" 


4 


0.458736 


0.389932 


0.304588 


8.09 


A233 = 


-4.06542 x 10" 


-5 


A 231 


= 7.40962 x 10" 


4 


0.459528 


0.388981 


0.304805 


8.09 



Table 1: A selection of LNV parameters and charged lepton mixings at SPSla which fit the 
neutrino masses and mixings for a normal hierarchy. All the points shown have \ 2 < 10~ 3 . The 
parameter sets marked with f are ruled out by the /xTi— > eTi constraints in [27]. 

As mentioned in section [|, the numerical analysis is performed in a basis where all 
quark mixings reside in the down quark sector. In the case of up type quark mixings, only 
A' couplings with the same quark flavour indices can be used to balance the contribution of 
A in the dynamical generation of the bilinear parameters, because the down quark Yukawa 
matrix Yrj is diagonal in this basis. Apart from this, the main effect in changing the quark 
mixing is expected to be a rescaling of A', with only minor changes to A and the charged 
lepton mixing angles. This is confirmed by our numerical results. 

In Tables || and|I], the branching ratios of the neutralino LSP for non-zero A233 and A233 
with the normal and inverted hierarchy respectively are generated using ISAJET7.64 [64]. 
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Ai 






A 2 






8 l 12 


#13 


#23 


Apt 


x 2 


a \f 
A 233 ~~ 


-5.69116 x 10" 


-0 


A233 


= 1.36023 x 10- 


4 


1.55779 


0.815115 


0.146126 


755 


0.01 


a \f 
A 233 — 


-5.69126 x 10" 


-6 


A211 


= 1.32365 x 10- 


-4 


1.38843 


0.760045 


0.140903 


758 


0.06 


a \ 1 
A 233 — 


-5.67940 x 10" 


-6 


A 123 


= 1.42938 x 10" 


-4 


1.81386 


-0.757538 


0.141975 


726 


0.05 


b \' — 
A 122 — 


-1.96283 x 10" 


-4 


A 122 


= 1.24824 x 10" 


4 


0.134765 


0.101938 


0.798282 


988 


0.43 


b \' — 
A 122 — 


-1.93673 x 10" 


-4 


A132 


= 1.47364 x 10' 


-4 


3.03234 


0.0866645 


0.931616 


743 


2.85 


b\ f — 
A 122 — 


-1.96175 x 10" 


-1 


A 123 


= 1.45708 x 10" 


-4 


0.144386 


0.0943266 


0.689708 


736 


0.52 


^311 = 


-1.90386 x 10" 


-3 


A311 


= 1.16446 x 10- 


4 


1.76438 


0.833990 


1.40752 


767 


1.41 


%n = 


-1.90027 x 10" 


-3 


A321 


= 1.44392 x 10- 


-4 


1.52379 


0.854066 


1.42182 


719 


0.22 


f A 3 ii = 


-1.90809 x 10" 


-3 


A 322 


= 1.37574 x 10- 


-4 


1.57829 


0.820421 


1.42557 


984 


0.01 


^233 = 


-1.36245 x 10" 


-4 


A 211 


= 1.90793 x 10- 


3 


1.55054 


0.815304 


0.146304 


979 


0.01 


^233 = 


-1.36235 x 10" 


-4 


\' 

A 221 


= 9.50337 x 10- 


-4 


1.55302 


0.815293 


0.146092 


975 


0.00 


^233 = 


-1.36016 x 10" 


-4 


\' 

A 231 


= 1.68685 x 10- 


-3 


1.56378 


0.815077 


0.145634 


753 


0.00 



Table 2: A selection of LNV parameters and charged lepton mixings at SPSla which fit the 
neutrino masses and mixings for an inverted hierarchy. The parameter sets marked with f are ruled 
out by the /xTi — » eTi constraints in [27]. Note that because of the high fine tuning, all digits 
present in the table are required to obtain the correct numerical results. 



In both cases there are many competing channels for the LSP decay, not surprising due to 
the typically large charged lepton mixing. On the other hand, it is amusing to observe the 
approximate left handed \i — r symmetry for the normal hierarchy, and approximate /U — r 
symmetry for both chiralities in the inverted hierarchy. This corresponds to an approximate 
symmetry among the A couplings after rotating to the diagonal charged lepton basis, and 
should be compared with the approximate fi — r symmetry in the experimentally favored 
near tri-bi maximal neutrino mass matrix [65]. This effect appears in all parameter sets 
we have checked which give rise to the observed neutrino mass pattern. However the result 
depends on the assumed condition of Zn, = Zir. We have checked that in other basis, for 
example with Zm = 13x3, this approximate symmetry ceases to appear. 





Channel 


BR 


Channel 


BR 


A122 




0.006 




0.006 


A132 


e~ v T n + 


0.007 


T~ v e n + 


0.007 


A123 


e~ t+ 


0.029 




0.029 


A133 


e~ v T t + 


0.034 




0.034 


A231 


fi~ v T e + 


0.005 


t~ e+ 


0.005 


A232 




0.027 




0.028 


A233 




0.138 




0.140 



Table 3: Decay channels of the LSP Xi a t 
SPSla, with A233 and A 233 at Mx given in 
the top row of Table [l], and a normal hierar- 
chy. Each process has a charged conjugated 
counterpart which is not shown in the table. 
Only branching ratios > 0.005 are displayed. 





Channel 


BR 


Channel 


BR 


A122 
A132 


er a* + 
e _ v T /i + 


0.063 
0.056 




0.063 
0.057 


A123 
A133 


er Pfj, t+ 
e~ v T t + 


0.067 
0.059 


V e T+ 
T~ V e T+ 


0.067 
0.060 



Table 4: Decay channels of the LSP \i a t 
SPSla, with A233 and A 233 at Mx shown in 
the top row of Table |^, and an inverted hier- 
archy. Each process has a charge-conjugated 
counterpart which is not shown in the table. 
Only branching ratios > 0.005 are displayed. 



Another observation is that the dominant decay channels in these two samples are 
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purely leptonic despite the presence of the A' couplings. This is mainly due to the choice of 
the third family indices in A233, which involves bottom Yukawa couplings in the dynamical 
generation of /ij. The cancellation condition then requires this coupling to be relatively 
small. In cases with A' coupling which involves smaller quark Yukawa couplings, the 3 body 
decay via this coupling can become dominant. Then, LSP decays into a charged lepton 
plus jets and some SUSY production events may be fully reconstructed since there would 
be no intrinsic missing energy. 



5. Discussion 



We have shown that it is possible to obtain the observed neutrino mass pattern with 2 
trilinear lepton number violating parameters in the mSUGRA set up in a basis where 
the lepton Yukawa coupling matrix is non-diagonal. However, depending on the flavour 
structure, the two non-zero couplings can either have one dominate over another, or have 
comparable magnitude. Here, we do not address the origin of hierarchies in the GUT-scale 
dimensionless couplings. They may be set, for example, by a spontaneously broken U(l) 
family gauge symmetry [66]. As an illustration, the first three parameter sets shown, la- 
belled with superscripts (a), in both Tables [l] and || have A233 a factor ~ O(10) smaller 
than the other LNV operator and so may be described by flavour models with only one 
dominant operator. For the next three parameter sets labelled with superscripts (b), Xuj 
and A' 12 2 are of similar magnitude. A rough indicator of whether particular LNV parame- 
ters may contribute significantly to the neutrino mass may be defined by noting that the 
dynamical generation of m scales like A(Y^) + 3X'(Y^) from Eq. ([2.17 ) . As an example 



suppose we are able to identify one coupling, say Xiku that contribute to the neutrino mass 
matrix significantly. Other couplings will also be important if 

AjfcZ _ (Yj5)mn /g 



Xjmn ( Y E)kl ' 

where the Yukawa couplings are evaluated in the weak interaction basis. Whether we 
have single coupling dominance or not can then be estimated by the ratio of the Yukawa 
couplings. However, for an inverted hierarchy, a slight change in the neutrino mass matrix 
can modify the results significantly, so more care is needed when applying this rough 
indicator. 

In the following we briefly compare our models with the literature. In neutrino models 
with bilinear R-parity violation, e.g. [17,18,30], the vacuum oscillation solution is generally 
favoured. This is due to the alignment effect which creates a large hierarchy between the 
tree level mass and the loop induced masses. Because we are interested in models with 
only trilinear R-parity violation at the unification scale, we make no further comments on 
this and refer interested readers to the literature. 

In [23], models with many A^- fe couplings are discussed. It was shown how such mod- 
els naturally lead to the now excluded vacuum oscillation solution of the solar neutrino 
problem. Our work complement such results by discussing the issues need to be addressed 
for a near tri-bi maximal mixing solution, and also performs the numerical calculation by 
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considering the full set of RGEs and include all one-loop contributions to the neutrino 
masses. 

In [33], fits on the neutrino masses in LNV mSUGRA using trilinear couplings in the 
diagonal lepton basis were attempted. The need for tree level mass suppression and the need 
to break the alignment effect were also discussed there. A hierarchy was assumed which 
included only the terms A^ 33 and/or Aj33 as a first approximation. This requires a different 
tree mass suppression method involving the RPC soft SUSY breaking parameters, whereas 
in our work, the suppression is obtained through interplay among the LNV parameters. 
The RGEs used in our work also included the contributions of the quarks and leptons from 
the first and second generations to the soft SUSY breaking terms. These can be important 
as, depending on their size, the LNV couplings with first and second generation indices 
may contribute significantly to the neutrino masses through coupling to the RPC Yukawa 
matrices in the RG equations. This also allows us to go beyond the natural hierachy 
assumed in [33]. Preference towards the small mixing angle MSW solution was preferred 
in their analysis, but only the mass squared differences and the atmospheric mixing angle 
was fitted. In our work, we obtain the best fit values of the LNV parameters with all 3 
PMNS mixing angles as well as the mass squared differences. 

In [25] , the authors described how the tri-bi maximal mixing pattern may be obtained 
by setting particular combinations of LNV parameters at the electroweak scale, which 
essentially involves setting two set of parameters, each with a different scale, to fit the 
observed neutrino mass pattern with the two scales lAm^l and Am^. Besides working in 
different basis, our models are in general inequivalent, because by setting the LNV couplings 
at Mz, the problem associated with a tree-loop hierarchy is absent, and the mass scales 
involved can be set directly, without any need of tree-loop level cancellation. Depending 
on the parameters chosen, it is even possible to have a vanishing tree level contribution in 
such models. The models presented in this paper are more difficult to control in this sense 
because the two mass scales involved are related by the (suppressed) dynamical generation 
of the bilinear terms. In a similar spirit, Ref. [67] assumes ansatzes for the weak-scale 
tri-linear LNV couplings in terms of six parameters. Some resulting models were shown 
to fit low-energy bounds as well as the atmospheric and solar neutrino anomalies. Some 
implications for LHC signals were investigated. 

Finally, we comment on the amount of tuning required. It may appear that the cancel- 
lation required is unnatural. To some extent we do need certain degree of tuning. We must 
cancel the tree-level mass sufficiently such that the loop corrections are roughly O(0. 1) 
times the tree level mass for a normal hierarchy. This implies that /ij is suppressed by 
0(1) compared to its 'natural' scale. This intuition is in agreement with the fine tuning 
measure discussed in section f|. For an inverted hierarchy, the loop corrections need to be 
roughly of the same order of magnitude as the tree level piece, and the solar neutrino mass 
difference arises from the difference between the masses of the two heaviest neutrinos, both 
of which are of the atmospheric mass scale. This means that while the parameters do need 
to be fine tuned for an inverted hierarchy due to the small solar neutrino mass difference 
compared to the atmospheric mass difference, for the normal hierarchy the tuning is rea- 
sonably mild. Once the tree and the loop contributions to the neutrino masses become 
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comparable, we see that the suppression of tree level mass also leads to departure from 
alignment, and significant neutrino mixing results if one allows for charged lepton mixings, 
possibly provided by some underlying flavour model. Thus large PMNS mixing is not 
unreasonable, although we have to treat the experimentally favoured near tri-bi maximal 
neutrino mass pattern as accidental in our approach. 

To summarise, we have investigated the minimal configuration in mSUGRA with tri- 
linear R-parity violation that can reproduce the observed neutrino oscillation observables. 
We obtain this by adding lepton number violating parameters in a basis where the charged 
lepton Yukawa matrix is not diagonal. We discuss on general grounds that it is impossible 
to obtain the neutrino mass hierarchy with just one LNV parameter in such a basis. This 
can be traced to the fact that the bilinear lepton number violating operators responsible 
for the tree level mass generation are radiatively generated by the trilinear coupling. De- 
spite the radiative origin of the bilinear operators, it turns out the tree level mass scale is 
dominant over the loop mass scale. Another general result is the alignment effect, whereby 
the tree level mass matrix and the loop corrections have almost the same pattern, therefore 
suppressing the light neutrino mass further. 

These effects can already be alleviated by introducing a second trilinear operator at 
Mx- The reason is that the bilinear operators will be suppressed if the trilinear operators 
contribute to their dynamical generations in an opposite way. The introduction of another 
dominant coupling also weakens the alignment, and the PMNS mixings observed can be 
obtained with large charged lepton mixings. 

We present, for both hierarchies, several data sets at SPSla obtained by minimizing 
the chi-square of the neutrino oscillation observables. The couplings involved are too small 
to affect the spectrum and decay of the sparticles, except for the LSP, which has to decay 
via the LNV operators. We find typically that there are a number of significant decay 
channels, while the actual channels are model dependent. Depending on the choice of the 
A' coupling, certain parameter sets are ruled out by the fi — e conversion in nuclei, but the 
LNV contributions to the branching ratios for other FCNCs considered are typically not 
experimentally accessible. 

While we need a certain degree of tuning to ensure the mass ratios fall in the observed 
range, the degree of tuning required for a normal mass hierarchy is relatively mild with 
Aft ~ O(10). The fine tuning required for the inverted hierarchy Apr ~ 0(800) reflects 
the need to arrange two a priori unrelated parameters to give two (quasi-) degenerate 
masses to account for the small solar neutrino mass squared difference. 

Evidently with so many lepton number violating operators, there are many ways to 
reconstruct the neutrino mass matrix. While we have shown that it is possible to build 
simple models with a manageable LNV parameter space, future work is required to realize 
the many other possibilities within this class of models. 
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A. Mass Matrices 

In this section we write down the CP-even, CP-odd, neutral fermion and charged fermion 
mass matrices. The notation follows that in [28], where the mass matrices are presented in 
terms of couplings in the Lagrangian. However following SDFTSUSY, the Standard Model 
Higgs vev is normalized to v ~ 246GeV. The 't Hooft Feynman gauge is used for numerical 
computation, although we present the mass matrices in the general £ gauge. 

Our convention is the same as in [25] where all mass matrices are presented in 2- 
component Weyl notation (apart from the inclusion of the sneutrino vevs). We also follow 
their definition of the mixing matrices. In the R-Parity conserving case, the mass matrices 
reduce to those in [41] 8 . 

With LNV, the down type higgs and the lepton doublets are put on the same footing. 
It is therefore instructive to write down the superpotential and the Lagrangian in a way 
that express this equivalence. We define 

C a = (H d ,Li), (A.l) 

where a, (3,... G {0, 1, 2, 3} and i, j, k, . . . G {1, 2, 3}. The parameters are defined by 

A Q /3fc = (^Ojk = (XE)jk, Kjk) , h a /3k = (flQjk = (tlE)jk,hijk\ 

Kjk = ( X 'ojk = ( Y D)jk,K jk ^J, h' ajk = (h' Qjk = (h D )jk,h' ijk J, (A.2) 



and 



= (( m c)oo = ml d , {m 2 c ) 0i = m\ iL .,{m 2 c )ij = (mf)ij). (A.3) 
In this notation, the superpotential is written as 

W = ^a/3k^a^i3Ek + \' a j k £ a QjDk + (Ytj)ijQiH u Uj - fi a C a H u , (A. 4) 

while the soft supersymmetry breaking Lagrangian is denoted to be 

—£soft = -^h a pk£a£pEk + h' a j k £. a QjDk — b a £ a h u + {hu)ijQih u Uj + h.c. 

+C\m\)C + E\ml)E + m\h\h u + Q f (m|)Q + W(ml)t) + D^ml)b 



+ 



^MiBB + ^M 2 WW + ^M 3 gg + h.c. 



(A.5) 



The mass matrices presented below follow the above notation. 

8 In Ref. [41] the mass matrices are presented in 4-component Dirac notation. 
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Higgs-sneutrino masses 

After electroweak symmetry breaking, the sneutrinos, v\ mixes with the Higgs bosons h u 
and h d = i> , resulting in CP-even (CPE) and CP-odd (CPO) scalars. The CPE and CPO 
Higgs-sneutrino mass eigenstates are obtained in a generic basis with D a = (/i^, v-i) after 
the diagonalization of the 5x5 mass matrices 



^ 'ylmh u Imi> 7 ^ M. 



2 

CPO 



( lmh u \ 
y Imi> 5 J ' 



where 



M 2 



CPE 



Vu 



+ 



g 2 + g 2 



M 2 



b a v a 



CPO 



i 9 2 + gh 



K -V U V5 \ 
-VuV^ VjVS J 

vl -v u v 5 

-V u Vj VjVS 



where £ is a gauge parameter 9 and 
Here we have \v a \ 2 = v\ + vf. The rotation matrices are defined by 



{g 2 + gl) ( 2 



(A.6) 
(A.7) 

(A.8) 
(A.9) 

(A.10) 



(A.ll) 
(A.12) 



z r m cpe z R = diag[m^ ,m^o,m?i], i = 1,2,3, 
Z a M cpo Za = diag[m G o,m A o,m-i ], i = 1, 2, 3. 
Neutrino-neutralino masses 

In a general flavour basis, with vp = (h? d , z^), the 7x7 neutrino-neutralino mass matrix is 
given by 

/ -iB \ 

-iW 3 



- (- 



iB - iW 6 K 



o 



M 



N 



(A.13) 



where (at tree level) 



M 



N 



\gv u - 

-\g2v u \g2V(3 



/Mi 

M 2 

\gv u \g2v u o -fip 



(A.14) 



The rotation matrix is given by 

Zj^M N Z N = diag[x°, n. 



a = 1,2, 3, 4. i = 1,2,3. 



(A.15) 



9 £ = 1 corresponds to the 't Hooft Feynman gauge. 



- 26 - 



Charged lepton-chargino masses 

The charged lepton-chargino mixing is displayed in the 5x5 matrix given below 

•iW~ e La )M c h+ \+h.c, (A.16) 




where ei a = (h d , e L ), and at tree level 

( M 2 -j=g 2 v u k \ 
Mc = K ^ ! . . (A.17) 

The diagonalisation matrices are defined by 

ZIM C Z+ = di ag [xt,ei], a = 1,2. i = 1,2,3. (A.18) 

B. One loop self energies of the neutral fermions 

For completeness, we list in this appendix the general form of the one loop contributions 
to the neutral fermion. They can be separated into 4 classes, depending on whether the 
chirality flip is in the loop or on the external legs, and whether the boson in the loop is 
a scalar or a vector boson. The complete listing in the RPC limit including the Feynman 
rules is in [41]. The Feynman rules for the general LNV calculation can be found in [25]. 
Our expressions are equivalent to the latter, but for the gauge boson contributions, we 
use the Passarino-Veltman (PV) 'B' functions [44] only, and present in a form easier to be 
compared with [41], which adopts the 't Hooft-Feynman gauge (£ = 1). 

Recall in the main text that the one loop self energy correction 5A4ij takes the form: 

SMij = {Y> D )i 3 - {M N )i k (Z L ) kj , (B.l) 

where Sd and Xj, are the mass and kinetic term corrections respectively. 



^ , \ ^0 ^0 (M N ) ik 




(a) contribution (b) contribution 

Figure 5: One loop self energies: scalar contributions. The et's and b's represent the vertex 
Feynman rules that can be found in [25,41]. 



The contributions from the scalar loops in Fig. [5] are given by 

167T 2 (fc 2 ) = -a^b^m^B (k 2 ,ml,ml), (B.2) 

167r 2 (££) fej (fc 2 ) = bl b^B 1 {k\m%,ml). (B.3) 
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Xi 



( M N)ik 



dt 



x° k i> 




(a) contribution 



(b) contribution 



Figure 6: One loop self energies: vector boson contributions. 




(B.4) 



(B.5) 



In the above expressions, the PV functions are defined as in [41]. It should be noted 
that the definitions of B\ in [41] and [25] differ by a sign. 

C. Mass insertion approximation in CPE-CPO cancellations 

It is well known that loop contributions from CP even (CPE) and CP odd (CPO) scalars 
tend to cancel each other in many physical processes, because the CPO scalars couple to 
other fields in the same way as the CPE scalars, apart from the extra i's associated with 
the Feynman rules. In particular, the radiative corrections to the neutrino masses involve 
strong cancellations among these contributions and may lead to numerical instability. The 
numerical fluctuation turns out to be irrelevant for the physically interesting LNV param- 
eter space we work in, but as we need to access regions where the light neutrino mass is 
typically suppressed, this issue needs to be addressed in order to obtain numerically stable 
results. Another situation where this may be important is when one wishes to consider 
heavy neutral scalar masses. In this case the numerical fluctuations may be comparable to 
the physical neutrino masses. 

Instead of working directly with the cancellation, we note that in the R-parity con- 
serving limit the cancellation between the CPE and CPO loops is exact. Therefore it is 
instructive to compute the deviation from this exact cancellation arises from the LNV con- 
tributions. We first concentrate on the CPE part. Eq. ( |Q| ) can be diagonalized in two 
stages, denoted by 




(C.l) 
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In the RPC limit, Z MIA will be a unit matrix. After the first stage of diagonalization, we 
get 



M 2 _^ I {M 2 H ) a 5 ab (Z°?* H Z9) aj , 

MCPE 1 (Zr a r Z H)ib h (C2) 



where {i, j} = 1 — 3, {a, b} = 1 — 2, and the hatted quantities are diagonal matrices. This 
matrix is then diagonalized by Z MIA , which can be obtained by successive approximations 
commonly known as mass insertion approximation (MIA). One advantage of performing 
the diagonalization in these two steps is that in the 't Hooft-Feynman gauge, the 3x3 
sneutrino part of the CPE and CPO matrices are the same, despite the fact that they also 
include even powers of lepton-number violating contributions. This means that the terms 
which contribute to the splittings between the CPE-CPO cancellations are in (ZjfanZp). 
This procedure also allows us to work in any flavour basis, in which the off diagonal entries 
M~ may be large. 

Now we can use MIA to obtain the mixing matrix Z MIA and the modifications to the 
eigenvalues. In general, for a symmetric matrix of the form 

M = M a 5 af s + e a ^ a) , (C.3) 

where e a p <C M a are matrix perturbations, the mixing matrix Z defined by 

Z T MZ = M (C.4) 

can be approximated to second order in e/M to be 

Z M = 8^ + (SZJfr + (SZ^fr + 0(e 3 /M 3 ), (C.5) 

where 

{SZ e ) Pl = e / 3 1 /(M J — Mp), 

{5Z e 2) Pl = ~E P e/3 p e n /2(Mp - M p )(M 7 - M p )5 Pl (C.6) 
+ E P W/3)/ 2 ( M 7 " M/j)(M 7 - M p ), 

and M is a diagonal matrix. From here it is easy to obtain an approximate correction to 
the eigenvalues 

SM a = M a - M a = Y -jgg^- ■ (C.7) 

This can be used to obtain the analytic expansion of Eq. (3.5) in [25]. In the basis in 
which we perform our diagonalization of M' N , this is given by 

s=l r=l 

\zR(2+i)s Z R(2+j)sBo(0, m 2 H o , m 2 K o) - Z A ^+i)s Z A(2+j) s Bo(0,'m 2 A o,m 2 .o) (C.8) 

where rn H o r and m A o r are the CPE and CPO neutral scalars, m K o ^ are the neutral 
fermions propagating in the loop, and Zr and Z A are the CPE and CPO scalar mixing 
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matrices respectively. Specifically, we seek an expansion of the expression in square brackets 
in Eq. (|Cj) 

Z R(2+i)s Z R{2+j) s Bo(0, m\o,m 2 K o) - z A(2+i)s Z A{2+j) s Bo(0, m \o,m 2 K o) 
~ Z R(2+i)h Z R{2+j)hBo(0,m 2 h0 ,m 2 K o) + Z R(2+i)H Z R{2+j)HBo(0,m 2 HO ,m 2 K o) 

- Z A(2+i)G Z A(2+j)G B o(0,m 2 :;0 ,m 2 o) - Z A{2+i)A Z A(2+j)ABo(0,m 2 AO ,m 2 K o) ^ 
+ ( Z R(2+i)v s Z R(2+j)v s - z A(2+i)u s z A(2+j)u s )Bo(0,rnlcp ,m 2 K o) 

+ (m 2 -cp E - ^lcpo) Z R(2+i)u s Z R(2+j)u a Co(^lcpo,m 2 -cPo,^l j Q) 
+ higher order terms. 

One can then substitute in the results from the MIA to obtain an analytic approximation 
to 0(a 2 /M 2 ). It is rather lengthy to show the full analytic expansion here. We merely 
note that in the limit of zero sneutrino vevs and no sneutrino mixing, the above expression 
agrees with Eq. (4.10) in [25], apart from the off-diagonal contributions of 5Z e 2, which is 
not included in Ref. [25] but is important for our purpose. 

D. Calculation of I 1 — > l J ^ 

In this section we give an explicit calculation of the branching ratio of leptonic FCNC: 
BRil 1 -> / J 7). For the phase space calculation, 4 component spinor convention as in [68] 
is used. For the calculation of the loop integrals, 2 component spinor notation is used 
instead. See the Appendix of [25] and references therein for details on the 2 spinor notation 
calculations. 

The effective Lagrangian for decay I 1 — > l J j is given by 

C eff = el J a^(A L P L + ArPr^F^, (D.l) 

where a^ v = §[7^,7^], and Pl/r = 5(1 T 75)- ^-L(R) are obtained from finite part of one 
loop diagrams to be described later. In this notation, the matrix element for the process 
is given by 

M = 2eu l ju^{K L P L + ArPr) Ui i€^ u , (D.2) 

where , Ujj are 4-component Dirac spinors. Summing over final state spins in the matrix 
element squared, and averaging over the initial I spin, we obtain 

k = - y \m\ 2 

2 z — / 

= 2e 2 Tr{a^(A L P L + ArPr)^ + m,/)(AJP fl + A* R P L ) 

v pa (rfiJ + mij)(-riw)Qvqo-} 
= 16e 2 (\A L \ 2 + \AR\ 2 yq- PlI )(q- Pl j). (D.3) 

The decay width of I 1 —> l J ~/ is given by the standard formula 

= l-mfr(\A L \ 2 + \AR\ 2 ), (D.4) 
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Using the result T(l — > l J u J i> 1 ) = G^m|/192-7r 3 (see for example Ref. [69]), the branching 
ratio is then related to BR^ 1 — > l J v J v I ) by 

BRil 1 _> = 1^!J(|A L | 2 + lAnfiBRil 1 -> Z W), (D.5) 

where the branching ratios BR^ 1 — > l J v J v I ) are given by [38] 

£i?(r -» /ii^V T ) = 0.1736 ± 0.0005 
Si?(r -> ez7 e z/ T ) = 0.1784 ± 0.0005 

BiJ^i -» ev e v^) ~ 1. (D.6) 

The calculation of the Wilson coefficients and can be performed in the 4- or 
the 2-component formalisms. To compute the loop integrals it probably is slightly easier 
to use 4 component conventions. On the other hand, the 2 component formalism is more 
transparent, especially to visualize the helicity flips from the fermion mass insertions. We 
will use the latter method. 



r J i 1 ■ 



1> 



mj mi 
l cJ l 1 --*-- 



icJ 



Figure 7: I 1 — > Z J 7 decay 

The loop integrals for are displayed in Fig. ^. For a loop with scalar (j> and fermion 
if), the couplings correspond to a left-handed, negatively charged lepton I 1 from the 
SU(2) doublet which flows into the vertex, whereas the bj^ couplings correspond to left- 
handed, positively charged anii-lepton l cI from the SU(2) singlet, again flowing into the 
vertex. The A# diagrams are similar to the Al ones, and can be obtained by interchanging 
a with b* and by reversing the helicity flow. Further details can be found in [25] and [56]. 

The analytic expressions for I — > l Jr y are given by: 
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where Q^, and are the electric charges of the scalar and fermion respectively, flowing 
away from the first vertex in the decay. n c is the number of colour in the loop. The loop 
integrals are given by: 

F li x ) = — ui 2 + 3x - 6x 2 + x 3 + 6xlnx], 
12(1 — xp 

F i( x ) = Tn7T~ 71 i 1 ~ 6x + 3x2 + 2x3 ~ 6z 2 hix], 

12(1 — xp 

Fs(x) = 1 [-3 + Ax - x 2 - 2lnx], 
2(1 — x) 6 

FAx) = s J l - x 2 + 2xlnx]. (D.9) 

2(1 — x) d 

Strictly speaking, the above loop integrals are only valid if the external fermion masses 
are negligible. Contributions from LNV terms typically involve SM fermions in the loop, 
so the external fermion masses cannot be neglected for a full calculation. However, because 
the loop scalar remains massive, it is possible to obtain corrections in powers of mf/m 2 ^, 
which are negligible for our purpose 10 . To be more explicit, we use the integral F% as an 
example. The full integral is given by 

i _ r 1 . r 1 -* . i 



771 



9 i<3 = / (lz / dy— 7) q ; r 7)— ; r tt . (D.10) 

' ' z(m^-yrnp-(l-y-z)mfj) + (l-z)m^ 



Thus the approximation in Eq. ( |D.9|) is valid up to ~ 0(mf /mV). Readers who are in- 
terested in applications of the full loop integrals, for example in a supersymmetric theory 
where the loop scalar mass can be light, may consult [70] for instance. 

We have checked explicitly that the above expressions are in agreement with [68] 
and [58] (in the context of b — > S7) 11 . We also checked that the expression for (g — 2)^ as 
well as the numerical results of the SPS benchmarks in the RPC limit agree with [71]. 
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